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1.0  INTRODUCTION 

This  is  the  final  report  prepared  by  Decision-Science  Applications, 
Inc.,  (DSA)  under  contract  number  F30602-78-C-01 29 , Space  Radar  Clutter 
Characterization,  sponsored  by  Rome  Air  Development  Center.  The  contents 
of  this  report  comprise  a model  for  the  description  of  radar  clutter  as  it 
is  involved  in  space-based  radars  (SBR).  It  is  expected  to  be  of  value  in 
modeling,  analysis,  and  design  of  future  SBR  systems. 

With  this  purpose  in  mind,  the  clutter  model  was  designed  to  be  as 
simple  and  easy  to  use  as  possible  while  providing  an  adequate  description 
of  clutter  as  it  is  expected  to  be  encountered  in  SBR.  The  present  data 
base,  based  on  measured  ground  clutter,  was  organized  into  seven  generic 
terrain  types.  Variation  of  the  mean  backscatter  coefficient  was  obtained 
by  regression  analysis  as  a function  frequency  and  grazing  angle  for  each 
terrain  type.  In  addition,  the  statistics  of  variation  about  the  mean 
were  explored.  Experimental  data  was  used  to  infer  a probability  density 
function  (pdf)  (sometimes  referred  to  as  a distribution)  and  a correlation 
length  for  several  of  the  terrain  types.  Using  a theory  of  scattering  for 
a number  of  independent  scatterers,  we  have  developed  a method  for  deriving 
the  pdf  as  a function  of  resolution  cell  size.  The  method  was  used  on  five 
terrain  types  to  determine  threshold  settings  to  give  a specified  false 
alarm  rate. 

This  report  is  organized  into  six  sections.  Following  this  intro- 
duction is  a section  which  discusses  the  general  problem  of  clutter  for 
look-down  radar,  pointing  out  the  special  considerations  for  the  SBR  case. 
Immediately  following  are  sections  which  present. in  order:  parametric 
regression  analysis  of  o° , statistics  of  variation  from  the  mean,  clutter 
spectra,  and  finally  recoirmendations  for  further  measurement  and  analysis. 
Also  included  is  an  appendix  describing  the  DSA  computer  program  which 
simulates  clutter  signals  based  on  a microscopic  description  of  clutter. 

This  was  developed  for  the  purpose  of  evaluating  signal  processing  algorithms 
for  realistic  clutter  scenarios  and  was  partially  supported  under  this 
contract. 
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2.0  GENERAL  DISCUSSION  OF  CLUTTER 

Before  examining  the  quantitative  details  of  the  model,  a qualitative 
overview  is  in  order.  This  section  describes  the  attributes  of  clutter 
with  emphasis  on  its  likely  impact  on  SBR  systems.  Figure  1 shows  a typical 
look-down  scenario,  graphically  illustrating  how  reflections  from  back- 
ground clutter  interfere  with  targets.  Within  thti  radar  footprint,  there 
is  distributed  clutter  such  as  vegetation,  discretes  such  as  buildings  and 
vehicles,  and  airborne  clutter  like  birds,  weather,  and  chaff.  Some  are 
stationary,  while  others  may  be  moving. 

In  the  missions  which  must  detect  air  vehicles  against  the  earth's 
background,  the  space  radar  is  confronted  with  the  traditional  look-down 
radar  problem  of  suppressing  ground  clutter.  The  space  radar  problem  is 
quite  different  from  the  airborne  radar  because  of  the  broad  extent  of 
clutter  in  the  range/range-rate  domain,  larger  resolution  cells,  and  range/ 
Doppler  ambiguous  returns. 

2.1  MAJOR  CLUTTER  SOURCES  AND  THEIR  IMPACT  ON  SYSTEM  PERFORMANCE 

Before  one  can  accurately  assess  the  adequacy  of  the  clutter  suppres- 
sion approach,  one  must  be  able  to  quantify  the  various  attributes  of  the 
clutter  background.  This  section  presents  a short  review  of  radar  clutter 
sources.  Those  aspects  of  clutter  which  most  impact  system  design  are 
presented  in  order  of  importance. 

The  basic  clutter  cross  section  determines  the  amount  of  clutter 
suppression  necessary  through  the  following  equation. 

i-mfs.  (2-D 

where  I is  the  required  AMTI  improvement  factor  required,  S/C  is  the 
necessary  signal -to-cl utter  ratio,  is  the  clutter  radar  cross  section 
(RCS),  and  is  the  target  RCS. 
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SATELLITE 


STRIP  OF  RANGE-RESOLVED 
BACKGROUND  CLUTTER 


Figure  1.  Schematic  of  Spaceborne  Look-Down  Radar 
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For  land  and  sea  clutter  aQ  is  determined  by  the  product  of  the 
backscatter  coefficient,  a0,  and  the  area  of  the  illuminated  resolution 
cell,  taking  into  account  the  additional  clutter  from  range  and  Doppler 
ambiguities.  In  general,  a0  is  dependent  upon  terrain  type,  grazing  angle, 
frequency,  and  polarization.  In  addition,  it  is  sensitive  to  such  factors 
as  moisture  content,  season,  time  of  day,  etc.  Sea  clutter  is  very  sensi- 
tive to  wind  speed  and  direction  and  wave  height  as  well.  Figure  2 shows 
some  typical  values  of  o°  as  a function  of  terrain  type.  Most  SBR  clutter 
is  of  the  extended  type,  however,  some  discrete  scatterers  will  be  present. 
They  lie  entirely  within  a resolution  cell;  hence  their  radar  cross  section 
is  sufficient  to  describe  their  impact. 

For  volume  distributed  clutter  sources  such  as  aurora,  weather,  birds 
and  chaff,  the  parameter  of  interest  is  the  RCS  density  per  unit  volume  n. 
The  clutter  RCS  is  given  by  the  product  of  n and  the  resolved  clutter  cell 
volume.  As  a general  rule,  these  sources  are  much  weaker  than  land  clutter 
but  can  appear  at  finite  velocities  with  respect  to  the  ground.  Figure  3 
summarizes  the  magnitudes  and  velocities  of  common  clutter  sources. 

2.2  GROUND  CLUTTER  DESCRIPTION 

In  addition  to  the  expected  value  of  clutter  RCS,  the  probability  dis- 
tribution about  the  mean  value  is  important.  The  probability  distribution 
impacts  the  required  S/C.  For  most  clutter  sources,  the  distribution  is 
similar  to  a Rayleigh  function.  However,  land  clutter  at  low  grazing  angles 
and  high-resolution  sea  clutter  has  been  shown  experimentally  to  exhibit 
a lognormal  distribution  of  large  variance.  The  impact  of  lognormal  clutter 
on  threshold  setting  for  detection  is  shown  in  Figure  4.  It  is  speculated 
that  the  distribution  approaches  Rayleigh  as  the  resolution  cell  size 
increases  for  all  types  of  clutter.  This  has  been  demonstrated  to  a certain 
extent  by  SKYLAB  X-band  dataJ 

Soofi,  Clutter  Model  for  Land,  Forest,  Snow,  Seaice,  and  Ocean,  Remote 
Sensing  Laboratory,  University  of  Kansas  Technical  Dept.  No.  TR-2923-2, 

July  1978. 
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Clutter  Distribution 


Of  interest  also  is  the  correlation  of  clutter.  This  spatial  corre- 
lation is  a measure  of  the  distance  between  independent  areas  of  clutter. 
Spatial  correlation  is  important  in  evaluating  the  tracking  performance 
and  in  design  of  constant-false-alarm  rate  detectors.  Temporal  correlations 
or  spectra  are  another  way  of  expressing  the  velocity  spread  of  clutter 
important  for  moving  target  indication  (MTI).  For  SBR  applications,  this 
is  usually  dominated  by  the  spread  induced  by  platform  motion.  However, 
the  inherent  internal  mot 'on  can  be  important  in  some  cases. 

2.3  CLUTTER  AS  VIEWED  BY  SPACE-BASED  RADAR 

Those  aspects  of  clutter  which  are  peculiar  to  space-based  radar  are 
summarized  below.  They  require  special  care  in  modeling  since  the  data 
base  for  SBR  clutter  is  limited. 

2.3.1  Large  Resolution  Cells 

Due  to  the  long  ranges  involved,  typical  footprints  measure  ten  to 
twenty  miles  by  as  much  as  one  hundred  miles.  The  geometry  is  such  that 
this  results  in  strips  which  extend  the  width  of  the  footprint  in  cross 
range  and  approximately  C/2BW  in  range,  where  C is  the  speed  of  light  and 
BW  is  the  modulation  bandwidth.  When  a pulse-burst  waveform  is  used,  there 
are  ambiguities  in  range-Doppler  space.  When  Doppler  ambiguities  are  spaced 
wide  enough  apart  to  prevent  blind  speeds,  range  ambiguities  occur  in  the 
mainbeam  footprint.  These  range  ambiguities  cause  clutter  foldover  in 
range,  adding  to  the  effective  clutter  per  resolution  cell.  These  ambigui- 
ties are  independent  in  amplitude  and  phase,  a fact  which  reflects  itself 
in  the  resulting  pdf. 

2.3.2  Discretes 

Because  of  the  large  resolution  cells  encountered  in  SBR,  it  is  almost 
inconceivable  (although  theoretically  possible),  for  a single  discrete 
scatterer  to  add  a significant  contribution  to  the  whole.  We  believe, 
since  our  model  is  based  on  data  from  actual  terrains  (not  set  up  in  a 
laboratory),  that  the  proper  contribution  of  discrete  scatterers  is  taken 
into  consideration.  This  is  true  for  parametric  regression  analysis, 
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variability  statistics,  and  spatial  correlations.  For  these  reasons  we  saw 
no  need  to  focus  explicitly  on  discretes. 

This  is  not  to  say  that  discretes  will  never  be  a problem  for  SBR. 

If  they  are  moving  at  velocities  which  overlap  the  range  of  possible  target 
velocities  and  have  an  RCS  comparable  to  targets  of  interest,  they  will 
look  like  targets.  These  may  have  to  be  acquisitioned  and  tracked  and 
eventually  rejected  based  on  a higher  order  detection  logic.  Evaluation 
of  the  impact  of  objects  of  this  sort  is  straightforward  but  not  well  defined 
enough  to  be  within  the  scope  of  this  project. 

2.3.3  Sidelobe  Clutter 

Sidelobe  clutter  is  a major  problem  for  SBR  as  it  is  for  airborne 

radars.  In  both  cases,  the  earth  is  illuminated  by  sidelobes  of  the  antenna 

pattern  as  well  as  the  mainlobe.  For  airborne  radars,  parts  of  the  earth 

illuminated  by  sidelobes  can  be  at  ranges  much  shorter  than  that  of  the 

a. 

target  and  mainlobe  clutter.  The  R'  term  in  the  radar  equation  can  signi- 
ficantly increase  sidelobe  clutter.  The  broad  velocity  spectrum  of  the 
sidelobe  clutter  compounds  the  difficulty. 

However,  for  SBR,  the  range  differences  are  not  so  great.  DSA  has 
the  capability  to  evaluate  sidelobe  clutter  parametrically  or  for  a specific 
antenna  pattern.  An  example  of  parametric  results  is  shown  in  Figure  5. 

A good  rule  of  thumb  is  that  sidelobe  levels  80  dB  (two-way  power)  below 
the  peak  are  sufficient. for  sidelobe  clutter  suppression  for  bomber-type 
targets.  Also,  sidelobe  levels  required  for  desirable  antijamming  proper- 
ties are  usually  much  lower  than  this.  For  smaller  targets,  the  sidelobe 
clutter  can  become  a more  severe  problem  than  mainlobe  clutter.  Accurate 
estimates  of  sidelobe  clutter  can  be  derived  from  the. model  if  they  become 
the  dominant  clutter  source. 
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Figure  5.  Effective  Clutter  for  -80  dB  Sidelobe  Level  as 
Function  of  Number  of  Doppler  Bins. 


3.0  MODEL  FOR  EXTENDED  CLUTTER 


3.1  GENERAL 

This  section  deals  with  a parametric  description  of  the  backscattering 
properties  of  extended  clutter.  In  particular,  the  backscatter  coeffi- 
cient a0  is  modeled  as  analytic  expressions  whose  coefficients  are  deter- 
mined by  best  fits  to  measured  physical  data. 

The  goal  is  to  derive  a model  which  will  be  useful  for  describing 
clutter  as  it  affects  the  performance  and  design  of  spaceborne  radar 
systems.  Although  a secondary  result  of  this  effort  is  certain  to  be  a 
better  basic  understanding  of  the  scattering  properties  of  land  and  sea, 
that  was  not  the  purpose  c*  this  study.  Therefore,  wherever  possible,  a 
wide  variety  of  terrains  were  included  within  each  generic  terrain  type,  and 
data  were  included  from  a variety  of  sources.  The  alternative  was  to  choose 
a single  experiment  or  source  upon  which  to  base  the  model.  We  feel  that 
the  former  approach,  although  it  may  not  provide  a tightly  controlled  data 
base,  provides  a better  representation  of  the  radar  community's  collective 
understanding  of  clutter.  This  also  was  the  approach  selected  by  I IT  in 
the  Overland  Radar  Technology  program.1 

An  essential  requirement  is  that  the  resulting  model  be  simple.  That 
is,  in  order  to  be  useful,  it  is  necessary  that  as  little  as  possible  need 
be  specified  a priori.  The  essence  of  the  physical  data  must  be  distilled 
into  a few  equations.  With  this  in  mind,  the  world  was  broken  down  into 
seven  terrain  types: 

1 . Desert 

2.  Rural  (fields,  orchards,  etc.) 

3.  Forest  (forest,  woods,  heavy  vegetation) 

4.  Mountains 

5.  Snow  and  ice 

o.  Urban  (cities  and  urban  areas) 

7.  Sea 

J.  Greens tei n,  et  al . , A Comprehensive  Ground  Clutter  Model  for 
Airborne  Radars,  I IT  Research  Institute,  Final  Report  Under  Contract  No. 
F33615-69-C-1387. 
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There  appeared  to  be  little  if  any  consistent  variation  in  o°  between 
Horizontal  and  vertical  polarization  for  land  clutter,  and  for  this  reason 
polarization  was  not  included  as  a model  parameter.  Thus,  grazing  angle 
and  frequency  were  the  only  parameters  other  than  terrain  type.  For  the 
case  of  sea  clutter,  fits  were  also  made  to  cross  polarized  data  (HV  and 
VH).  Also,  sea  data  were  separated  into  categories  on  the  basis  of  wind 
speed.  Specifically  they  were  separated  into  two  groups,  less  than  or  equal 
to  15  kt  and  greater  than  15  kt,  and  fit  separately. 

The  fitting  was  made  easier  by  developing  a simple  computerized  data 
entry  and  data  base  retrieval  system.  The  data  and  the  system  will  be 
retained  after  the  completion  of  this  project. 

3.2  RESULTS 

The  process  of  arriving  at  an  analytical  model  for  distributed  clutter 
involved  the  following  steps: 

1.  Select  and  organize  a data  base. 

2.  Preview  the  data  to  isolate  the  most  important  parameters. 

3.  Select  analytic  forms  which  appear  to  be  likely  candidates 
for  providing  good  fits. 

4.  Select  the  best  fits  by  experimentation. 

5.  Eliminate  data  which  deviates  too  m>;ch  from  the  norm. 

6.  Perform  the  final  fits  to  determine  the  adjustable  coefficients. 

7.  Repeat  steps  1 through  6 if  necessary. 

The  analytic  form  which  was  most  successful  in  fitting  the  land  data 
was  as  follows: 

5°  = A + BQ  + (C  + De)  log  F (3-1) 

where  we  have  adopted  the  notation  5°  which  is  a0  expressed  in  decibels  (dB) 
(5°  = 10  log-^  a0),  and  where  e is  the  grazing  angle  in  degrees  and  F is  the 
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frequency  in  GHz.  Forms  involving  terms  proportional  to  F,  log  e,  and  sine, 
were  tried  and  found  to  be  less  suitable.  For  sea  clutter,  on  the  other 
hand,  5°  falls  off  sharply  at  low  grazing  angles  and  it  was  discovered  that 
the  grazing  angle  dependence  was  better  expressed  by  a term  proportional 
to  log  e.  Hence 


5°  = A + B(log  e)  + C(log  F)  + D('!og  e log  F)  (3-2) 

The  method  of  least  squares  or  regression  analysis  is  well  suited  to 
these  expressions.  The  result  is  an  analytic  expression  5°(e,  F)  such 
that  the  sum  of  the  squares  of  the  error  terms  for  N data  points 


S = 


I CS’<VFf)  - s?]2 

i~l 


(3-3) 


is  minimized.  This  of  course  does  not  imply  that  the  sum  of  the  squares 
of  the  error  terms  is  minimized  if  a°  is  not  expressed  in  decibels.  In 
fact,  one  might  argue  that  a least  squares  fit  to  5°  is  always  biased 
slightly  low  for  this  reason.  Therefore,  we  compute  for  each  fit  an  adjust- 
ment factor,  which  we  recommend  be  added. 

5°  = o°  + 6 
2 

where  6 = 0 . 1 1 5 ( S/N ) where  <5  is  given  by  Equa.  3-3.  In  short,  this  cor- 
rection is  arrived  at  by  assuming  that  the  variation  of  the  data  about  the 
fit  is  Gaussian  (in  dB).  The  correction  term,  5,  is  the  difference  in  dB 
between  the  mean  o°  and  the  mean  5°. 

The  remaining  details  of  the  results  and  the  procedures  through  which 
they  are  arrived  are  discussed  below,  seoarately  by  terrain  type. 
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3.2.1  Desert 


The  types  of  terrain  included  under  the  broad  category  of  "desert" 
ranged  from  "sandy  desert"  to  "dry  lake  bed."  The  data  as  a whole  were 
very  well  described  by  the  fit  below,  the  square  root  of  the  mean  squared 
error  being  3.2  dB. 

a°  = -34.15  + 0.1630  + (8.69  + 0.00936)  log  F (3-4) 


with 


6 = 1.16  dB 

Figures  6a  and  6b  show  plots  of  this  formula  versus  e and  versus  F, 
with  and  without  data  superimposed. 

3.2.2  Rural 

Of  all  the  data  gathered  on  the  radar  reflection  from  land  clutter 
the  vast  majority  is  based  on  crops  and  farmland.  Some  data  included 
initially  had  to  be  excluded  because  it  deviated  too  much  from  the  norm 
and  adversely  affected  the  fits.  Examples  of  excluded  data  were  bare 
plowed  fields  and  Arizona  farms.  In  the  end,  a fit  based  on  data  from 
five  experiments  yielded  a fit  which  was  a fair  description  of  rural 
clutter  data.  This  was 

5°  = -23.61  + 0.9946  + (3.53  + 0.0916)  log  F 
6 = 0.79  (3-5) 

Plots  are  shown  in  Fig.  7a, b.  The  root  mean  square  (rms)  deviation  from 
the  fit  was  2.62  dB. 

3.2.3  Forest  (Heavy  Vegetation) 

Unlike  the  previous  case,  there  is  relatively  little  data  available 
on  forest  clutter.  Nevertheless,  data  was  found  sufficient  to  cover  the 
frequencies  and  grazing  angles  of  interest. 
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Figure  6a.  Desert  Terrain,  d°  versus  Grazing  Angle 
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Figure  6b.  Desert  Terrain,  a6  versus  Frequenc. 


Much  of  the  data  used  in  the  fit  comes  from  the  University  of  Kansas 
Remote  Sensor  Laboratory  (RSL),1  who  performed  a separate  curve  to  fit 
their  data.  Their  fit  does  not  predict  so  low  a 5°  at  low  frequencies 
and  low  grazing  angles  (-22.9  dB  at  e = 10°,  F = 0.5  GHz).  However,  their 
fit  was  linear  in  F,  not  log  F and  it  noticeably  deviates  from  the  observed 
data  at  low  frequency  and  low  grazing  angles. 

The  chosen  analytical  form  fits  the  data  quite  well  with  an  rms  differ- 
ence = 2.81  dB.  The  result  was 

5°  = -42.36  + 0.520  + (24.93  - 0.358e)  log  F 
6 = 0.91  dB  (3-6) 

This  is  plotted  in  Fig.  8a, b.  Note  that  in  the  region  of  low  grazing 
angles  and  low  frequencies  a very  low  5°  is  predicted.  At  frequencies 
below  1 GHz  there  was  no  data  for  grazing  angles  below  60°.  Therefore, 
this  result  should  be  treated  with  caution. 

3.2.4  Mountains 

Understandably,  there  is  not  a great  deal  of  mountain  data  available. 

Also  there  is  a great  variability  in  the  radar  reflection  from  mountains 
due  to  the  variability  in  the  slope  and  composition  of  the  terrain,  making 
parametric  analysis  difficult.  Mountains,  nevertheless,  are  of  great 
interest  to  SBR  analysis  due  to  the  possibility  of  large  specular  returns 
from  facets  normal  to  the  incident  radiation. 

Most  of  the  data  were  of  Arizona  mountains  taken  by  Naval  Research 
Laboratory  (NRL).^  There  were  a few  measurements  taken  by  the  Environ- 
mental Institute  of  Michigan  (ERIM),^  from  L-band  calibrated  SAR  data. 

Again  the  data  were  fairly  true  to  form  and  the  rms  deviation  from  the  fit 
was  2.98  dB.  The  fitted  parameters  are 

R.  K.  Soofi,  ibid. 

^J.  C.  Daley,  NRL  Terrain  Clutter  Study,  Phase  II,  NRL  Report  No.  6749, 

Oct  1968. 

A.  Maffatt  et  al.,  L-Band  Radar  Clutter  Statistics  for  Terrain  and  Ice, 

ERIM  Report  No.  128900-9-F,  (I),  January  1978. 
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GRAZING  ANGLE  <DEG> 


FREQUENCY  <GHZ> 


= -30.24  + 0.1736  + (8.82  - 0.0156)  log  F 


(3-7) 


6 = 1.02  dB 

This  is  plotted  in  Fig.  9a, b.  The  data  show  that  5°  continues  to 
decrease  as  e decreases.  However,  due  to  the  specular  scattering  phenomenon 
and  shadowing,  data  from  mountains  and  hills  tend  to  have  a very  large 
variance.  This  could  mean  that  the  clutter  could  actually  be  more  severe 
a problem  as  the  grazing  angle  decreases. 

3.2.5  Snow  and  Ice 

Snow  and  ice  covered  terrain  is  of  great  interest  in  the  planning  and 
analysis  of  SBR  early  warning  systems  since  the  northern  perimeter  surveil- 
lance fence  covers  arctic  regions.  There  exists  a fairly  large  amount  of 
data  on  ice  and  snow. 

There  was  a high  degree  of  variance  in  the  data  from  the  different 
sources  so  the  rms  deviation  from  the  fit  is  larger  than  in  the  previous 
examples,  being  5.01  dB.  The  fitted  parameters  were 

5°  = -32.97  + 0.3409  - (1.797  + 0.035e)  log  F 
6 = 2.9  dB 

with  3°  plotted  in  Fig.  10a, b. 

One  curious  feature  of  this  data  is  that  5°  decreases  with  frequency. 
Again,  Kansas,  Remote  Sensoring  Laboratory*  performed  independent  fits  to 
their  snow  data  covering  grazing  angles  from  20°  to  90°  and  1 to  18  GHz.  In 
general,  the  frequency  dependence  was  weak  in  their  fits;  for  some  categories, 
5°  decreased  with  frequency  and  for  others  5°  increased. 


^R.  K.  Soofi,  et  a!.,  ibid. 


i 

\ 


'A'*'  — 


gure  9a.  Mountain,  5°  versus  Grazing  Angle 


GRAZING  ANGLE  <DEG> 


Figure  10b.  Snow  and  Ice,  00  versus  Frequenc, 


3.2.6  Urban  Terrain 


Urban  terrain,  which  covers  cities  and  residential  areas,  differs 

from  the  other  types  of  terrain  in  that  the  reflected  energy  is  probably 

dominated  by  that  from  man-made  objects,  for  this  reason,  there  is  a great 

deal  of  variability  of  the  data  from  a given  source  as  well  as  from  one 

source  to  the  next.  For  example,  data  on  the  city  of  Phoenix J show  a much 

higher  5°  at  0 = 8°  for  all  frequencies  than  any  other  grazing  angle.  This 

was  excluded  from  the  fit  after  it  appeared  to  be  adversely  affecting  the 

2 

e dependence.  Data  on  the  city  of  Chicago  were  excluded  because  the 
frequency  dependence  was  not  monotonic  and  varied  greatly. 

In  the  end,  a fit  was  obtained  with  an  rms  deviation  of  3.8  dB.  The 
parameters  were  as  follows 

5°  , -12.42  - 0.0600  - (1.79  - 2.07e)  log  F (3-9) 

6 = 1.65  dB 

This  is  plotted  in  Fig.  11a, b. 

3.2.7  Sea  Clutter 

Compared  to  land  clutter,  sea  clutter  is  relatively  benign.  It  is 
very  important  that  a model  for  sea  clutter  be  included  in  our  SBR  clutter 
model,  however,  because  most  of  the  earth  is  covered  by  ocean  and  one  would 
like  to  know  how  small  a target  can  be  detected  in  the  sea  clutter  environ- 
ment. 


Due  to  the  importance  oi  sea  clutter  to  navy  radar  operations  and  to 
theoretical  interest,  there  exists  a good  deal  of  data.  There  is  a fairly 
strong  dependence  on  such  factors  as  wave  height  and  wind  velocity  all  of 


^J.  C.  Daley,  et  a!.,  fb’d. 

2W.  Ament,  "Radar  Terrain  Reflections  for  Several  Polarizations  and 
Frequencies,"  Trans.  1959  Symp.  Radar  Return,  Pt  2.,  May  1959. 


31 


Urban,  o°  versus  Grazi 


which  are  not  independent.  There  is  also  a stronger  polarization  dependence 
than  in  the  case  with  land  clutter.  For  these  reasons,  we  have  segregated 
that  data  into  four  categories  based  on  wind  speed  and  polarization. 

These  are  as  follows: 


Category 

Polarization 

Wind  Speed 
(kts) 

1 

like 

>15 

2 

like 

<15 

3 

cross 

>15 

4 

cross 

<15 

For  a single  sea  clutter  model,  category  one  is  the  preferable  choice 
since  like  polarization  is  more  common  and  high  winds  are  prevalent  in 
northern  seas.  Sea  clutter  tends  to  drop  sharply  at  low  grazing  angles  and 
for  this  reason  better  fits  resulted  when  the  e dependence  was  proportional 
to  log  e instead  of  e.  The  results  of  the  fits  are  given  below  for  each 
category. 


Category  1 . (like  polarization/high  wind) 

a3  = -58.9  + 19.5  log  9 + 1.19  lgo  F + 3.60  log  e log  F 
ms  deviation  = 5.3  dB 

6 = 3.2  dB  (3-10) 

Category  2.  (like  polarization/low  wind) 

a°  = -65.5  + 21.96  log  0 + 6.68  log  F = 1.45  log  9 log  F 
rms  deviation  = 6.9  dB 

6 = 5.5  dB  (3-11) 
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L (cross  polarization/high  wind) 


5°  = -59.8  + 12.3  log  0 +1.43  log  F + 0.88  log  F log  0 
rms  deviation  = 1.5  dB 


6 = 0.25  dB 


(3-12) 


(cross  polarization/low  wind) 


5°  = -58.8  + 9.62  log  0 + 5.68  log  F - 7.69  log  0 log  F 
rms  deviation  = 2.2  dB 
6 = 0.56  dB 


(3-13) 


These  are  plotted  respectively  in  Figs.  12a  through  15b. 

Summary 

The  fitted  parameters  for  the  model  as  well  as  the  rms  deviations 
and  6,  the  factor  which  compensates  for  fitting  j°  rather  than  g°  are 
presented  in  Table  1.  In  general  we  feel  that  these  fits  provide  a fair 
description  of  the  collective  clutter  data  available  at  this  time. 


GRAZING  ANGLE  <DEG> 


Figure  12b.  Sea.  Like  polarization,  high  wind  speed,  5°  versus  frequency. 


, low  wind  speed,  a°  versus  frequency. 
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Sea.  Cross  polarization,  high  wind  speed,  a°  versus  grazing  angle 


TABLE  1 

FITTED  PARAMETERS  IN  REGRESSION  ANALYSIS 


TERRAIN  TYPES 

A 

B 

C 

D 

<*  (dB) 

7 

rms 

deviation 

(dB) 

Desert 

-34.15 

0.163 

8.69 

0.0093 

1.16 

3.20 

Rural 

-23.61 

0.0994 

3.53 

0.091 

0.79 

2.62 

Forest 

-42.36 

0.520 

24.93 

-0.358 

0.91 

2.81 

Mountains 

-30.24 

0.173 

8.82 

-0.015 

1.02 

2.98 

Snow  and  Ice 

-32.97 

0.340 

-1.797 

-0.035 

2.9 

5.01 

Urban 

-12.42 

-0.060 

-1.79 

2.07 

1.65 

3.8 

Sea 

*1 

-58.9 

19.5 

1.19 

3.60 

3.2 

5.3 

2 

-65.5 

21.96 

6.68 

1.45 

5.5 

6.9 

3 

1 

cn 

LO 

00 

12.3 

1.43 

0.88 

0.25 

1.5 

4 

-58.8 

9.62 

5.68 

-7.69 

0.56 

2.2 

* 

Note  A,  B,  C,  D,  are  defined  differently  for  sea  clutter.  See  Eqs.  3-1 
and  3-2. 


4.0  PROBABILITY  DISTRIBUTION  OF  CLUTTER  o° 

Of  primary  importance  in  assessing  the  impact  of  clutter  on  detection 
is  not  the  average  clutter  cross  section  but  the  probability  that  it  will 
exceed  a critical  threshold.  This  notion  can  be  derived  from  the  probability 
density  function  (pdf)  for  RCS.  For  example,  most  theories  of  detection 
model  the  noise  contribution  to  the  in-phase  and  quadrature  components  as 
Gaussian.  This  implies  that  its  amplitude  is  Rayleigh  distributed 


P(V) 


(4-1) 


in  the  probability  that  the  amplitude  V is  between  V and  V + dV.  Likewise 
the  power  n is  exponentially  distributed 


P(n) 


(4-2) 


where 


Thus,  the  probability  that  noise  of  mean  NQ  will  exceed  n is  given  by 

-n/N 

Pc(n)  - e 0 (4-3) 

These  statistics  also  apply  in  theory  tc  the  reflections  from  a few 
ideal  clutter  types,  such  as 
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1. 


An  infinite  number  of  scatterers  in  each  resolution  cell 
(like  chaff) . 

2.  Any  number  of  scatterers  in  the  resolution  cell  if  their 
amplitude  distribution  is  Rayleigh. 

In  each  case,  the  phase  must  be  random  and  uniformly  distributed,  and  the 
amplitude  and  phase  of  each  scat.terer  must  be  independent. 

A Rayleigh  distribution  is  very  convenient  since  P,  the  probability 
density  function,  contains  only  one  parameter.  Then,  if  the  pdf  describing 
the  target  statistics  likewise  has  only  one  parameter,  the  mean  RCS  (e.g., 
nonfluctuating,  Swerling  I,  II,  III,  IV),  the  detection  statistics  are 
completely  determined  by  the  singal-to-interference  ratio. 

Unfortunately  for  the  detection  problem,  ground  and  sea  echos  do  not 
always  obey  Rayleigh  statistics,  thus  complicating  the  theory  and  making 
detection  in  general  more  difficult.  The  probability  distribution  of  clutter 
amplitude  has  received  a significant  amount  of  attention,  both  theoretically 
and  experimentally,  in  the  past.  Nevertheless,  there  is  still  a great  deal 
of  uncertainty  and  misunderstanding.  The  fact  is  that  variability  in  the 
clutter  reflectivity  is  very  difficult  to  define,  much  less  to  predict  or 
to  measure. 

In  this  section  we  shall  attempt  to  define  clutter  statistics  which 
are  relevant  to  the  SBR  detection  problem,  identify  factors  which  contribute 
to  clutter  variability,  and  show  how  one  might  use  the  model  to  predict 
the  probability  that  clutter  will  exceed  a given  threshold.  Throughout, 
the  size  of  the  resolution  cell  will  be  stressed  as  a crucial  parameter. 

In  the  case  of  SBR,  we  are  primarily  interested  in  large  resolution  cells. 

4.1  DEFINITION 

First,  we  shall  define  what  we  feel  to  be  the  problem  of  most  interest 
in  detection  in  clutter.  Generally,  one  is  interested  in  the  number  of 

*'  * i 
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false  alarms  per  unit  time  or  false-alarm  rate.  Given  the  number  of 
resolution  cells  per  second  this  implies  a false  alarm  probability.  Most 
modern  radars  contain  constant-false-alarm  rate,  CFAR,  capability  which  makes 
use  of  a local  average  to  set  the  detection  threshold  which  preserves  the 
false  alarm  rate.  The  CFAR  must  take  into  account  the  variability  in  the 
clutter  return  from  this  local  sample.  The  variability  of  clutter  in  a 
given  scene  is  important  in  determining  how  flexible  the  CFAR  must  be  in 
processing  a single  scene,  and  finally  the  variability  from  scene  to  scene 
is  also  important.  By  scene  we  mean  a batch  of  data  upon  which  detection 
decisions  are  made.  Typically,  a scene  would  include  all  the  data  within 
the  beam  footprint. 

4.2  CAUSES 

We  shall  list  a number  of  factors  which  influence  clutter  variability 
and  the  scales  over  which  we  feel  the  influence  is  most  prominent.  We  shall 
discuss  each  in  more  detail. 

4.2.1  Mutual  Interference  Effects 

Mutual  interference  effects  are  responsible  for  the  so-called  Rayleigh 
statistics  encountered  in  chaff  and  most  distributed  clutter.  The  pdf  is 
well  known  and  is  given  by  the  Rayleigh  theory 

, -o°h' 

P(o°)  = -f  e (4-4) 

where  p is  the  mean  value  of  o° . It  applies  only  to  cases  where  the  indi- 
vidual scattering  center  themselves  are  Rayleigh  or  for  an  infinite  number 
of  scatterers. 

4.2.2  Man-made  Objects 

We  can  use  urban  data  to  determine  the  cross  section  per  unit  area  of 
man-made  objects  by  assuming  the  urban  clutter  is  dominated  by  man-made 
objects.  The  difference  between  urban  clutter  and  rural  clutter,  for  example 
is  10  dB  at  low  grazing  angles  and  decreases  as  the  grazing  angle  increases 
(Section  3.2).  If  one  assumes  that  a fraction  p of  resolution  cells  contain 
all  man-made  objects  then  the  fraction  p will  be  10  dB  higher  than  the  others 
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As  resolution  cells  increase  in  size,  every  resolution  cell  will  contain 
part  natural  and  part  man-made  clutter,  and  the  variance  will  be  smaller. 
Taken  to  the  limit  of  very  large  resolution  cells,  perhaps  several  miles 
square,  one  might  expect  every  resolution  cell  to  contain  a fraction  p of' 
man-made  clutter.  Then  the  variance  should  reduce  to  essentially  zero. 

4.2.3  Variation  in  Terrain  Type 

In  a given  scene  (or  even  resolution  cell),  there  may  be  more  than  one 
terrain  type.  For  instance,  forest  and  fields  or  rural  and  urban  terrains 
often  coexist.  For  land  clutter,  the  difference  in  a0  among  terrain  types 
is  about  20  dB,  the  difference  between  urban  clutter  and  desert  or  snow. 
Again  the  actual  variance  is  a function  of  resolution  cell  size  with  the 
same  argument  as  given  in  4.2.2. 

4.2.4  Variation  ir  the  Slope  of  Terrain 

In  cases  where  the  terrain  is  not  flat  and  the  period  of  undulation 
is  large  compared  to  resolution  cell  dimensions,  the  actual  grazing  angle 
will  vary  from  cell  to  cell.  Since  clutter  reflectivity  with  few  excep- 
tions. is  strongly  dependent  on  grazing  angle,  this  leads  to  significant 
variance  in  the  observed  a°. 

This  hypothesis  is  easy  to  check.  If  we  approximate  the  terrain  by  a 
series  of  ridges  which  have  a sinusoidal  undulation,  Fig.  16,  the  probabil- 
ity density  function  of  grazing  angle  and  hence  RCS  can  be  computed  as  a 
function  of  the  ratio  of  amplitude  to  period  of  undulation.  A computer 
program,  which  included  terrain  masking  effects,  was  written  to  do  this. 
Qualitatively,  the  0.05  case  resembles  gentle  rolling  hills  and  the  0.1  case 
resembles  very  hilly  terrain. 

Figs.  17  and  18  show  the  pdf's  for  the  probability  tM'  the  actual 
local  grazing  angle  is  between  0 and  0 + d0  for  nominal  grazing  angles  of 
5,  10,  30,  50  and  90°  for  A/T  = 0.05  and  0.1,  respectively.  For  instance, 
in  Fig.  17  we  note  that  for  a nominal  grazing  angle  of  30°  and  A/T  equal  to 
0.05,  the  actual  grazing  angles  fall  between  11°  and  48°!  Furthermore, 
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ERRAIN  ROUGHNESS  MODELS 


Figure  16.  Model  for  Local  Terrain  Undulation 


0.05 


50 


Figure  17.  Probability  Density  Function  Versus  9 for  Undulating  T< 


0.05 


Figure  18.  Probability  Density  Function  Versus  0 for  Uridulati 


values  at  the  extremes  of  the  interval  rf  possible  values  occur  with  higher 
probability. 

Given  the  pdf's  for  grazing  angle  it  is  simple  to  convert  to  pdf's  for 
RCS.  This  is  done  in  Figs.  19  and  20  for  rural  clutter  at  3 GHz.  These 
figures  show  the  probability  per  dB  that  RCS  will  fall  between  5°  and 
5°  + do0.  Again,  for  the  30°  nominal  grazing  angle  o°  falls  between  -19.6  dB 
and  -14.4  dB  with  higher  probability  near  the  extremes.  Table  2 gives  the 
maximum  spread  in  c°  for  each  combination  of  nominal  grazing  angle  and  A/T 
ratio.  Note  that  at  low  grazing  angles  there  is  a lot  of  shadowing  taking 
place  so  that  there  is  a higher  probability  of  no  clutter  at  all.  All  of 
these  results  are  consistent  with  the  observed  phenomenon  that  the  ratio  of 
mean  to  median  clutter  increases  significantly  as  the  grazing  angle  decreases. 

4.2.5  Variation  in  Local  Conditions 

By  local  conditions,  we  mean  such  factors  as  moisture  content  of  soil, 
type  and  density  of  trees,  heights  of  buildings,  etc.  Perhaps  the  best 
estimate  of  this  variation  comes  from  the  variation  in  observed  measurements 
from  highly  controlled  experiments.^  From  our  data  base  we  have  have  esti- 
mated this  variation  for  each  terrain  type  in  Table  3. 

The  University  of  Kansas,  RSL,  has  collected  data  to  explicitly  show 
seasonal  diurnal  variations  in  clutter.  These  are  presented  in  Figs.  21  and 
22.  Other  data  showing  variation  of  5°  with  row  orientation  relative  to 
radar  show  large  variation  at  high  grazing  angles  (>50°)  but  the  variation 
appears  to  decrease  substantially  to  within  the  3 dB  estimate  given  in 
Table  3 at  lower  grazing  angles.  Figure  21  supports  a seasonal  variation, 
which  might  -150  be  interpreted  as  moisture  and  crop  maturity  variation,  of 
about  4 to  6 dB.  This  figure  also  supports  a variation  of  about  10  dB  with 
crop  type. 
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Soofi,  ibid. 
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Figure  19.  Probability  Density  Function  Versus  o°  for  Undulating  Terrain.  A/T  = 0.05 


TABLE  2 

SPREAD  IN  5° 

DUE  TO  UNDULATION 

OF  TERRAIN 

MAXIMUM  SPREAD 

IN  5° 

A/T 

5° 

2.3  dB 

0.05 

3.8  dB 

0.05 

10° 

5.2  dB 

0.05 

30° 

5.1  dB 

0.05 

50° 

2.6  dB 

0.05 

vo 

o 

o 

4.0  dB 

0.10 

5° 

4.0  dB 

0.10 

10° 

8.7  dB 

0.10 

30° 

9.5  dB 

0.10 

cn 

o 

o 

4.6  dB 

0.10 

90° 

TABLE  3 

ESTIMATED  5°  VARIATIONS 


TERRAIN  TYPE 

A o° 

Crops 

6 dB 

Forest 

5 dB 

Snow/' Ice 

2 dB 

Urban 

7 dB 

Sea 

3 dB 
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Figure  21.  Examples  of  Seasonal  Variation 

From  W.  Simkins,  private  correspondence 


Crop  Type:  Hybrid  of  Sorghum  (S  and  M) 
Polarization:  HH 
Frequency  (GHz):  7.25 


21  24 

Midnight 


3 6 9 

Time  of  Day 


12  15 

Noon 


Figure  22.  Examples  of  Diurnal  Variation 

From  W.  Simkins,  private  correspondence. 


57 


4.3  CLUTTER  AVERAGING  AND  RESOLUTION  CELL  SIZE 

It  is  an  established  fact  that  for  sea  clutter,  the  distribution 
differs  qualitatively  and  quantitatively  as  a function  of  resolution  ceil 
size.  In  particular,  if  the  dimensions  of  wave  height  and  period  are  com- 
parable or  smaller  than  the  resolution  dimensions,  the  distribution  tends 
to  a Rayleigh  formJ  There  is  experimental  evidence  to  support  the  same 
type  of  trend  for  land  clutter,  but  the  situation  for  land  is  complicated 
by  its  larger  inherent  variability  compared  to  sea  clutter. 

We  have  a measure  of  the  distribution  of  RCS  in  the  limit  of  extremely 

2 

large  resolution  cells.  This  is  provided  by  the  University  of  Kansas  RSL 

who  analyzed  SKYLAB  data  over  the  U.S.  and  portions  of  South  America.  The 

2 

resolution  cell  size  is  400  to  1000  km  . This  is  greater  than  the  correla- 
tion length  for  any  local  undulation  or  any  local  terrain  variation.  RSL 
did  observe  a finite  distribution  of  a°,  however,  even  after  the  inter- 
ference effects  were  eliminated  by  incoherent  averaging.  Two  examples  are 
shown  in  Figs.  23  a and  b. 

Kansas  explained  the  presence  of  long  tails  in  the  distribution  in 
terms  of  specular  scattering  from  large  lakes  and  Utah  salt  flats.  At  high 
grazing  angles,  a specular  return  results  in  a large  echo  producing  a tail 
at  high  5°.  At  low  grazing  angles,  the  specular  is  highly  forward,  result- 
ing in  a low  a°  for  backscattering  and,  hence,  a long  tail  at  the  lower 
side.  The  small  but  finite  width  of  the  central  peak  in  Fig.  23  is  due  to 
the  variability  of  the  terrain  from  cell  to  cell.  This  small  width,  in 
view  of  the  many  sources  of  variability  observed  for  smaller  cells,  is  cer- 
tain to  be  due  to  the  averaging  effect  of  the  large  cells.  It  is  interest- 
ing to  note  the  trend  in  SKYLAB  data  toward  broader  distributions  with 

^G.  V.  Trunk  and  S.  F.  George,  "Detection  of  Targets  in  Non-Gaussian  Sea 
Clutter,"  IEEE  Transaction  on  Aerospace  and  Electronic  Systems,  Vol . 

AES-6  No.  5,  September  1970. 

2S.  M.  Purduski,  "Distribution  Tests,' 5-Percentile  Values,  and  Autocorrela- 
tion Analysis  of  Skylab  S-193  Overland  Radar  Data,"  RSLTM  2923-1.  (1978) 
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GRAZING  ANGLE.  W'-SB0 
UPPER  DECILE-  -6.0  dB 
MEAN:  -7.33  dB 
MAXIMUM:  -5.0  dB 


POLARIZATION 
INCIDENCE  ANGLE 
DATA  SOURCE 
NUMBER  OF  SAMPLES 
MEAN 


V 

BETWEEN  32.0’  AND  34.0’ 
SL2  SOUTH  AMERICA 
4601 
-7.32 


Figure  23b.  Histogram  of  the  Scatterometer  Backscatter  Coefficient, 
South  America;  From  Sobti,  ibid. 
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decreasing  grazing  angle,  Fig.  24.  This  behavior  is  consistent  with  other 

'i 

data  also. 

The  averaging  effect  of  resolution  cell  size  is  complicated  by  other 
phenomena  which  depend  on  cell  size.  For  example,  some  terrain  may  resemble 
a geometrical  shape  with  a smooth  surface.  When  the  resolved  areas  are 
smaller  than  the  size  of  the  piece  of  terrain,  the  returns  are  from  indi- 
vidual pieces  and  may  appear  fairly  uniform.  However,  as  the  cell  size 
increases  the  piece  of  terrain  acts  as  a single  scattering  center  and  the 
distribution  of  observed  RCS  is  qualitatively  different.  An  example  of 
this  is  the  specular  parts  of  the  SKYLAB  distribution  discussed  above. 

Another  hypothetical  example  is  terrain  made  up  of  flat  areas  randomly 
oriented  with  respect  to  the  radar.  Each  facet  has  a radiation  pattern 
of  the  form  sin  (x)/x  where  x is  proportional  to  the  angle  of  the  radar 
relative  to  the  facet  nominal.  Assuming  x is  a random  variable,  the  distri- 
bution which  results  resembles  lognormal  when  each  facet  is  sampled  sepa- 
rately. 

i 

5 

5 

With  these  compl icatio’is  in  mind,  we  present  other  distributions. 

First  there  are  synthetic  aperture  radar  (SAR)  data1  measured  by  the 
Environmental  Research  Institue  of  Michigan  (ERIM).  The  SAR  technique  is 
used  normally  for  high-resolution  imaging  of  the  ground.  Figure  25  shows 
X-band  data  for  mountains,  crops,  and  a city,  for  two  resolution  sizes, 

10'  x 10'  and  20'  x 20'.  These  cells  were  obtained  by  incoherently  adding 
even  higher  resolved  cell  size. 

I 

2 

Another  example  involves  a RADC  study  of  the  Seek  Igloo  Radars.  ] 

These  are  ground  radars  whose  beams  can  be  directed  at  nearby  clutter.  f 

RADC  measured  probability  distributions  by  setting  thresholds  and  counting 

'ERIM,  private  communication. 

W.  L.  Simkins,  et  al . , Seek  Igloo  Radar  Study.  Rome  Air  Development  Center, 
RADC-TR-77-338,  (October,  1977),  A047897. 
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Figure  24.  Distribution  of  SKYLAB  Data  versus  Grazing  Angle 
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a0(dB) 

Grazing  Angle  15° 
(o  = 0.151  = -3.2  dB) 


Figure  25.  ERII1  X-band  SAR  Histograms  (Pol  = HH) 
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false  alarms.  Presumably,  Rayleigh- type  interference  effects  are  included. 
Additionally  there  is  a greal  deal  of  terrain  masking  due  to  the  low 
grazing  angles.  Resolution  cells  are  quite  large  and  were  varied  by  . 
changing  the  range  resolution  (or  pulse  width  x).  Several  terrain  types 
were  observed.  Most  were  hills  and  mountains  but  some  swamp  land  was  also 
included.  Examples  from  the  RADC  report  are  shown  in  Figs.  26  and  27.  In 
general,  the  distribution  from  hills  and  mountains  is  much  broader  than  the 
swamp  land.  Also  in  almost  every  case  the  probability  of  exceeding  a high 
value  of  5°  is  smaller  for  the  larger  resolution  cells  (e.g..  Fig.  26) 
indicating  that  an  averaging  is  taking  place.  The  extremely  broad  dis- 
tribution for  hilly  terrain  should  probably  be  replotted  with  the  terrain- 
masked  areas  eliminated.  This  is  shown  in  Fig.  28;  for  curve  20  in  Fig.  27, 
note  that  this  distribution  is  much  more  in  line  with  what  one  would  estimate 
based  on  the  variability  factors  we  have  discussed  above,  and  with  other  data. 

4.4  MANY  SCATTERER  THEORY 

We  have  discussed  clutter  distributions  so  far  in  a very  general  way, 
avoiding  until  now  an  analysis  of  its  impact  on  threshold  settings,  detection, 
end  false  alarms.  We  feel  that  too  much  has  been  made  of  the  shapes  of 
distributions  and  curve  fitting  to  predict  detection  statistics  and  false 
alarms.  The  reason  is  that,  for  acceptable  false  alarm  rates,  the  signal - 
to-cl utter  ratio  must  be  very  high;  therefore  only  the  extreme  tails  of  the 
distribution  are  of  real  interest.  There  is  no  data  at  this  extreme  so 
that  predictions  are  usually  based  on  excessive  extrapolation. 

Our  model  is  based  on  the  idea  that  the  expected  values  of  cr°  can 
vary  due  to  the  various  factors  listed  above,  and  that  one  can  estimate  and 
predict  a distribution  of  mean  a0. 

Let  us  define  amplitude  A as  So5".  The  problem  reduces  to  predicting 
amplitude  and  phase  of  the  combined  echos 


(4-5) 
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TERRAIN  TYPE  AREA  NO. 


REFLECTIVITY  (dB) 

(Tatalina),  drom  W.  Simkins  et  al 


CURVE  NO.  PLACE  DATE  t usee  TERRAIN  TYPE  AREA  NO. 

19  TATALINA  25  AUG  1 MOUNTAINS  9 


-40  -35  -30  -25  -20  -15  -10  -5  0 

REFLECTIVITY  (dB) 

Figure  27.  Comparative  Distributions  (Tatalina),  from  W.  Simkins,  et  al . 


Assuming  A.  and  $•  are  statistically  independent  and  $.  is  uniformly  distri- 
j J 1 .1  _ I 

buted  over  a full  period.  Beckmann  has  shown  that  to  order  N where  N 

is  the  number  of  independent  scatterers  the  pdf  of  V is  given  by 


P(V)  = 


2V 

I 


-V2/I 


1 + Jr 


£(A4) 

E2(A2) 


2V 


+ 1 


21* 


(4-6) 


where  I = NE(A^)  and  E(x)  denotes  expected  value.  Therefore,  defining 
2 0 

0 = V /IQ,  the  probability  of  exceeding  any  given  value  a is  given  by 


Pc(a)  = e (1  + 3[oV2  - a]) 


(4-7) 


where 


6 - (c  - 2) 


and 

a = E(a2)/E2(a)  = E(A4)/E2(A2) 

This  tells  us  exactly  how  the  averaging  process  of  large  cell  sizes 

influences  the  probability  distribution.  For  example,  assuming  a corre- 

? 

lation  length  i,  N approximately  equals  AR/x,  where  AR  is  the  resolution 
cell  area.  Eq.  4-6  shows  that  as  AR  approaches  infinity,  P(V)  approaches 
Rayleigh. 


In  the  limit,  where  resolution  cells  are  small,  N is  small  so  that 
each  cell  may  contain  a different  clutter  of  a different  expected  amplitude. 
In  the  case  of  N = 1 , then  obviously 

P(V)  = P(A) 


V.  Beckmann,  Probability  in  Communication  Engineering,  New  York:  Harcourt, 
1967. 
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Note  that  the  parameter  6 in  Eq.  4-7  is  actually  a measure  of  how 
closely  the  clutter  pdf  resembles  a Rayleigh  distribution.  For  6=0 
the  distribution  is  Rayleigh.  Also  note  in  the  definition  of  6 in  the 
line  below  Eq.  4-7  that  a small  value  of  6 can  occur  for  either  of  two 
reasons:  large  number  N of  independent  scatterers  in  the  resolution  cell, 
and  a = 2 which  occurs  if  the  underlying  pdf  of  individual  scatterers 
is  Rayleigh. 
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Figure  29.  Probability  of  Exceeding  Median  Versus 


Figure  30.  Plot  of  BETA  versus  N for  Lognormal  and  Log  Rectangular 


The  fact  is  that  there  is  no  single  number  such  as  correlation  length 
which  can  completely  describe  the  effect.  For  example,  there  is  a corre- 
lation length  associated  with  terrain  undulation,  another  associated  with 
changes  in  vegetation,  another  associated  with  major  variations  in  terrain 
type,  etc. 

The  Skylab  data^  shows  that  some  of  these  correlations  lengths  may 
be  very  large.  There,  the  resolution  cells  were  very  large,  400  to  1000 
square  kilometers.  Yet  there  is  a variation  of  about  five  dB  in  the  mean 
o°  from  cell  to  cell.  This  is  an  indication  of  either  the  variation  in 
scattering  properties  of  the  clutter,  or  a failure  to  average  enough  samples 
within  each  cell  to  remove  the  Rayleigh  effects. 

On  the  other  hand,  there  is  a component  which  has  a very  short  corre- 
lation length,  namely  discretes.  In  general,  one  could  probably  estimate 
the  most  important  correlation  lengths  by  looking  at  contour  maps  or  (even 
better),  SAR  radar  images. 

7 

Some  experimental  data  has  been  gathered  by  I IT  to  determine  the 
autocorrelation  function  of  clutter.  To  be  meaningful,  the  Rayleigh  effects 
should  have  been  averaged  out  and  the  resolution  of  the  experimental  radar 
should  be  less  than  the  correlations  of  interest. 

The  question  which  now  arises  is  how  to  estimate  the  number  of  inde- 

4 2? 

pendent  scatterers,  N,  and  P^(A)  from  which  to  compute  E(A  )/E  (A“). 
Logically  these  should  be  obtained  from  single-look  experimental  pdf's 
and  autocorrelation  functions  where  the  experimental  resolution  is  fine 
enough  to  resolve  individual  scatterers.  The  I IT  results  are  based  on  a 
100  ft2  resolution  area  and  some  of  the  Rayleigh  effects  are  averaged  out. 

If  P^( A)  is  based  on  many-look  samples,  and  does  not  include  the  Rayleigh 
spread,  one  can  define  P^,(A')  and  redefine  3'  as 

]S.  M.  Purduski,  ibid. 

2S.  Kayel  et  si . , Extensions  to  the  ORT  Clutter  Model,  I IT  Technical  Report 
under  Contract  No.  F33615-69-C-1387  11971). 
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6' 


(4-9) 


= _3_  £(A4) 

8N  E2(A2) 

If  not  all  of  the  Rayleigh  effects  ha'.e  been  averaged  out,  the  error  which 
results  from  redefining  0 in  this  way  makes  the  analysis  slightly  conserva- 
tive. Also,  in  general,  N1  < N since  averaging  tends  to  smooth  the  terrain 
return  resulting  in  a longer  correlation  length. 

The  effective  value  of  N‘  is  obtained  from  the  autocorrelation  func- 
tion r(x).  We  derive  N1  from  r in  the  following  manner: 

For  independent  samples,  the  ratio  of  the  variance  of  the  individual 

? 2 
samples,  a , to  samples  of  ensemble  averages  of  N samples  each,  4^,  is  just 

N.  Thus 


(4-10) 


o 

If  they  are  not  independent,  we  compute  4N  and  define  the  effective  Nas  in 
the  above  equation.  In  order  to  do  this,  it  is  convenient  to  confine  our- 
selves to  the  variable  component  of  the  clutter  return.  The  total  clutter 
return  is  given  by 


C(X)  = V(X)  + u 


where 


V(x)  dx  = 0 


(4-11) 


Therefore  the  following  discussion  is  simplified  by  considering  a function 
of  mean  zero,  v(x) . 
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If  V(x)  is  sampled  at  a large  number  of  x values  and  averaged  within 
an  interval  <Sx,  the  variance  of  the  averages  is  given  by 


2 

,1  _ 2iT 
' x<$ 


5 x 

r{x)  dx 


where  r(x)  is  defined  as  the  expected  value  of  the  product 


(4-12) 


V ( x 1 + x)V(x') 


divided  by  E[lr(x)]  and  is  computed  by  taking  the  autocorrelation  of  V(x) 


(4-13) 


I IT  computed  the  autocorrelation  functions  of  the  variable  component 
of  a°  for  four  terrain  types.  These  are  shown  in  Fig.  31.  These  are  only 
experimental  approximations  to  the  true  autocorrelation  functions  which 
should  be  smooth.  Our  approximations  to  the  true  values  based  on  I IT's 
experimental  evidence  is  shown  in  Table  4. 


Note  that  6x  is  a strip  of  continuous  clutter.  If  there  are  ambiguous 
range  strips  spaced  more  than  about  2 nmi  apart,  each  strip  contributes 
additional  independent  samples. 


Thus,  the  effective  values  of  N for  each  clutter  type  is  given  below 
from  Eq.  4-10,  for  Namb  ambiguous  ranges: 


N (effective) 

Desert  (<5x/0.79)  x Namb 
Farms  (Sx/0.72)  x Namb 
Mountains  (6x/0.24)  x Namb 
Cities  ( 6 x / 1 -12)  x Namb 
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Figure  31.  Autocorrelation  of  Variable  Component  of  Several  Terrains, 
from  S.  Kazel  et  al. 
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TABLE  4 

AUTOCORRELATION  FOR  FOUR  TERRAIN  TYPES 

TYPE  1.  r(x)  = a e"x/C  + b x < xQ 
= 0 x > xc, 

a b c 

( nmi ) 

Desert  0.8  0.2  0.18 

Farm  0.75  0.25  0.18 

Mountain  1.0  0.0  0.12 

1 

! 

| 

For  cities  the  form  is  slightly  different:  \ 

TYPE  2.  r(x)  = 0.6  e"X/0,27  + 0.4(1  - x/2)  x < 2 nmi  ! 

! 

=0  x > 2 nmi  > 


1.25 

0.9 


CO 


I IT  did  not  compute  r(x)  for  woods,  but  a correlation  length  of  1 nmi,  based 
on  other  terrains,  seems  a reasonable  and  conservative  estimate.  Therefore, 
we  estimate  for  woods, 

Neff  = (Sx/1)  X Namb 

4.6  EXAMPLES  AND  PREDICTIONS  FOR  FIVE  TERRAIN  TYPES 

This  section  describes  a procedure  for  predicting  the  threshold  set- 
tings, above  median  clutter,  required  to  achieve  a given  false  alarm  proba- 
bility. It  is  based  on  experimental  data  and  the  distribution  suggested  by 
Beckmann.  Examples  are  worked  out,  based  on  the  SAR  data  of  IItJ  This 
data  was  analyzed  expressly  for  the  purpose  of  clutter  modeling.  Five 
terrain  types  were  considered:  desert,  woods,  farms,  mountains  and  cities. 
Grazing  angles  were  not  reported. 

We  shall  begin  with  a review  of  the  theory  of  Sec.  4-4  in  order  to  set 
the  requirements  for  RCS  and  scattering.  It  states  that  if  amplitude  and 
phase  are  independent  and  the  pdf  of  amplitude,  P ( A)  is  independent  and 
uniformly  distributed,  then  the  pdf  of  the  power  resulting  from  the  sum  is 
given  by  Eq.  4-6  or  4-7  with  (3  replaced  by  ($'  of  Eq.  4-9  if  Rayleigh  effects 
have  been  eliminated  by  averaging. 

The  pdf's  for  desert,  woods,  farm,  mountain,  and  city,  are  shown  in 
Fig.  32.  For  desert,  woods,  and  mountain,  the  distributions  are  nearly 
lognormal.  For  farms  the  distribution  is  approximately  the  sum  of  two 
normal  distributions  or  contaminated  lognormal.  For  cities  it  is  well 
described  by  breaking  the  region  into  two  sections  and  performing  straight- 
line  fits  (on  log-log  paper)  in  each  sector.  The  resulting  parameters  are 
given  in  Table  5.  I IT  could  not  calibrate  the  data,  so  the  mean  values  of 
these  distributions  are  meaningless  in  themselves. 

^S.  Kazel  et  al . , ibid. 
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Figure  32.  Comparison  of  Probability  Density  Function:for  Various  Terrains 
from  S.  Kazel  et  al. 


TABLE  5 

PROBABILITY  DENSITY  FUNCTIONS 
1 and  2.  Desert  and  Woods 

P(u)  = — 2 e‘u2/2a2 

/2F  a 

u = 10  1cg10(A2) 
a = 2.1  dB 

3.  Farms 

P(u)  = — e-u2/2o2  + c-u2/202 

a /2T  a2 

o]  = 1.17  dB 
o2  ~ 8.60  dB 

4.  Mountain 

P(u)  = — 2 e“u2/2cj2 

/2tt  0 
a = 3.02  dB 


5.  City 

P(u)  = 3.4  x 10"15  u11,4 
= 9.23  u"1'9 


u < 15 
u > 15 


TABl.E  6 

SUMMARY  OF  STATISTICAL  PARAMETERS 


TERRAIN 

pdf 

E(A4)/E2(A2) 

N: 

e' 

Ti 

T2 

Desert 

lognormal 

1.27 

253 

0.0019 

12.6 

14.3 

Woods 

lognormal 

1.27 

200 

0.0024 

12.6 

14.3 

Mountain 

lognormal 

1.64 

833 

0.0022 

12.5 

14.3 

Farms 

contaminated 

lognormal 

1.69 

278 

0.0023 

12.6 

14.3 

Cities 

power  laws 

9.04 

178 

0.0190 

14.3 

15.5 

Assumptions: 


Sx 

^amb 


T 

T 


1 

2 


10  nmi 
20 

implies  10"®  false  alarm  probability 
-8 

implies  10"  false  alarm  probability 
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From  these  pdf's  we  compute  the  ratio  E(cr  )/E  {0),  and  hence  3 1 for 
each  terrain  type.  This  was  straightforward  for  all  terrain  except  "city". 
Here  the  second  moment  is  infinite  unless  an  upper  bound  is  placed  on  the 
highest  a.  We  chose  36  dB  because  this  is  the  highest  value  on  IIT's 
reported  pdf.  These  are  presented  in  Table  6 along  with  estimated  values  of 
effective  N'.  Also  using  Fig.  29  or  Eq.  4-7  we  obtain  the  appropriate  thres- 
hold, relative  to  the  median  a°,  required  for  the  given  alarm  prob«:  ilities. 

The  above  is  appropriate  for  a single  terrain  type  in  the  footprint. 

If  a CFAR  is  used  which  operates  on  local  areas  smaller  than  the  footprint, 
it  is  appropriate  for  the  more  general  case  where  each  local  area  is  com- 
posed of  at  single  type  of  terrain. 


We  have  considered  the  case  also  where  there  are  M areas  each  having 
a different  type  of  terrain.  Then,  there  is  roughly  a log  rectangular 
distribution  fcr  P(u)  and  the  number  of  independent  samples  is  N.  The  width 
of  the  distribution  is  a most  20  dB.  We  obtain 


E(A^i  = 
E?(A2) 


2.3 


Assum  ng  N ~10, 

3'  = 0.086 


We  conclude  from  this  exercise  that  there  is  little  cause  for  concern 
for  high  tail  clutter  returns  in  a typical  pulse-Doppler  surveillance  SBR. 

In  a SAR  mode,  however,  the  fine  resolution  would  greatly  reduce  the  number 
of  independent  scatterers  and  probably  require  thresholds  higher  by  1 to 
5 dB  than  those  called  for  by  Rayleigh  statistics.  Over  cities,  the  thres- 
hold would  have  to  be  10  to  20  db  higher  than  called  for  by  Rayleigh  statis- 
tics. 


Nevertheless,  although  the  pdf's  measured  by  I IT  cover  much  higher 
tails  than  most,  it  is  desirable  to  examine  the  high  tails  more  closely  to 
increase  our  confidence  in  3'. 
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5.0  CLUTTER  SPECTRA 

In  addition  to  the  total  magnitude  of  clutter  in  the  resolution  cell, 
the  frequency  spectra  of  its  radar  echos  is  also  critical  for  moving  target 


indication  (or  detection),  MTI.  This 
earth  reflection  will  totally  dominate 

the  Doppler  signature  must  be  used  to 
stationary  earth. 

However,  a number  of  factors  can 
Doppler  spread,  resulting  in  a cut-off 
These  factors  include,  platform  motion 
of  clutter,  and  glint  effects. 

Platform  motion  and  earth  rotati< 
shift  of  the  center  line  also  cause  a : 
The  spread  due  to  platform  is  given  by 


s because,  in  a look-down  geometry,  the 
most  tactical  targets.  As  a result, 

separate  moving  targets  from  the 

cause  the  clutter  to  have  a finite 
velocity  for  detecting  slow  targets, 
and  earth  rotation,  internal  motion 

n,  in  addition  to  producing  a finite 
pread  due  to  the  finite  beamwidth. 


P 8 cos  4 
x & 


(5-1) 


where  is  the  effective  beamwidth,  and  4 is  the  squint  angle. 

Earth  rotation  is  usually  a much  smaller  effect  but  can  b^rme  compar- 
able to  platform  motion  at  high  altitude.  In  general,  it  can  either  add  to 
or  cancel  the  spread  due  to  platform  motion.  In  a geosynchronous  orbit  they 
cancel  exactly. 

The  internal  motion  of  clutter  is  usually  very  small  for  wind-blown 
vegetation.  It  would  not  be  an  important  consideration  unless  platform 
motion  effects  were  eliminated  by  a displaced  phase  center  antenna  (DPCA) 
and  very  slow  target  detection  were  required. 

Other  types  of  clutter  with  internal  motion  are  weather  and  land 
vehicles.  Wind-blown  weather  has  a shift  relative  to  the  ground  and  a 


finite  width.  Land  vehicles  can  also  be  a problem  for  slow  target  detection 
in  that  they  have  an  RCS  of  the  same  order  of  magnitude  as  air-vehicle 
targets  and  there  may  be  many  of  them  present  in  each  resolution  cell.  Of 
course  this  is  less  of  a problem  in  the  arctic.  Table  7 gives  some  rough 
order  of  magnitude  velocity  spreads  from  various  sources.  Figure  33  is  an 
example  of  a measured  clutter  spectrum  from  wind-blown  vegetation. 
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Figure  33.  L-Band  Power  Spectral  Density  of  Alaskan  Land  Clutter 
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6 . 0 RECOMMENDATIONS  FOR  FURTHER  DEVELOPMENT 

We  feel  that  about  as  much  meaningful  information  as  possible  has 
been  extracted  from  the  present  clutter  data  base.  We  have  performed 
parametric  fits  to  seven  terrain  types  using  two  parameters  and  data 
from  several  sources.  In  addition,  the  University  of  Kansas  has  taken 
extensive  data  on  crops,  snow,  sea  ice,  and  forest.  They  have  performed 
detailed  fits  in  incidence  angle  and  frequency  for  different  polariza- 
tions and  time  of  year. 

However,  there  are  shortcomings  in  the  present  data  base  which 
could  be  improved  by  further  experiment  and/or  computer  synthesis.  The 
following  are  areas  where  more  information  is  needed: 

1.  arctic  terrains 

2.  low  grazing  angles  { < 10°) 

3.  large  resolution  cells 

4.  better  calibration  of  data  for  all  frequencies  and  terrains 

5.  distributions 

Data  is  needed  for  arctic  terrains  because  a surveillance  fence  is  required 
over  northern  Canada  and  Greenland  against  a Soviet  attack  against  CONUS. 
There  are  little  data  for  grazing  angles  less  than  10°  for  any  terrain 
type,  and  what  data  there  are  sometimes  contradict  each  other.  It  is 
dangerous  to  extrapolate  from  higher  angles  because  the  nature  of  the 
scattering  is  qualitatively  different  for  very  low  grazing  angles. 

The  space  radar  systems  popularly  considered  at  this  time  have  large 
resolution  cells.  That  is,  they  are  large  compared  to  the  resolution  of 
most  clutter  experiments  but  small  compared  to  that  of  SKYLAB.  To  test 
the  results  of  our  false  alarm  analysis,  resolution  cells  of  the  order 
of  10  miles  by  100  meters  with  a number  of  range  ambiguities  is  required. 
There  is  a need  for  better  calibration;  calibration  should  be  a major, 
factor  in  considering  future  measurements.  Finally,  more  emphasis  should 
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be  put  on  distributions  of  RCb  about  the  mean  values.  It  is  the  fraction 
of  clutter  cells  which  exceed  a threshold  several  dB  above  the  mean  which 
is  most  important  for  most  SBR  missions.  The  analysis  of  this  report 
indicates  that  the  distribution  is  effectively  Rayleigh  but  that  remains 
to  be  shown  experimentally.  Most  of  these  requirements  point  to  an  actual 
SBR  experiment;  however,  much  could  be  done  as  well  by  an  airborne  experiment. 

Three  candidate  experiments  have  been  considered: 

1.  Analysis  of  the  L-band  SAR  data  from  SEASAT. 

2.  Calibration  and  analysis  of  L-band  and  X-band  data 
taken  by  ERIM  over  various  sites  in  Canada. 

3.  An  airborne  experiment  designed  by  Spectra  Research 
Corporation  (SRC)  for  the  purpose  of  simulating  an 
SBR  clutter  scenario. 

Table  8 summarizes  these  three  experiments  with  respect  to  frequency, 
terrain,  grazing  angles,  and  resolution.  Of  the  three,  the  one  which  seems 
most  likely  to  produce  information  valuable  to  clutter  modeling  of  SBR  is 
the  ERIM  data.  With  SEASAT,  the  grazing  angles  are  too  high  and  there  is 
only  one  frequency.  Nevertheless,  SEASAT  data  does  exist  and  there  are 
plans  by  the  Defense  Advanced  Research  Projects  Agency  and  Lincoln  Laboratory 
to  analyze  it  with  emphasis  on  large-cell  clutter  and  targets.  Although 
the  SRC  experiment  covers  the  most  relevant  scenarios,  it  has  higher 
technical  and  financial  risk  because  the  required  hardware  has  not  been 
assembled. 


; '•  e is  also  the  possibility  that  computer  synthesis  of  data  will 
answer  some  of  the  pressing  questions  at  a fraction  of  the  cost  of  an 
actual  experiment.  This  would  be  accomplished  by  modeling  a scene  in 
great  detail  within  a variety  of  terrain,  hills,  discretes,  etc.  The 
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RCS  of  each  scatterer  would  be  estimated  from  extensive  existing  data.  Then 
the  pdf  of  clutter  magnitude  for  each  scene  could  be  generated  as  a function 
of  resolution  by  combining  the  scatterers,  in  amplitude  and  phase,  in  the 
computer  as  they  would  combine  in  a radar. 

These  are  a few  of  the  possibilities.  It  is  probably  safe  to  say  that 
clutter  will  never  be  understood  completely  because  there  are  so  many 
different  variables  to  control  (e.g.,  resolution,  frequency,  angles,  season, 
time  of  day,  etc.). 
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7.0  SUMMARY  AND  CONCLUSIONS 

In  summary,  a clutter  model  has  been  developed  for  use  in  space- 
based  radar  analysis.  Those  aspects  of  clutter  unique  to  or  prevalent  in 
space  radar  have  been  identified  and  discussed.  In  short,  those  aspects 
identified  were:  large  resolution  cells,  large  footprints,  range  ambiguous 
mainbeam,  high  platform  velocity,  and  (except  for  high  resolution  SAR) 
dominance  of  discretes  by  distributed  clutter. 

Existing  clutter  measurement  data  were  entered  into  a computerized 
data  base  and  categorized  by  terrain  type,  polarization,  frequency,  grazing 
angle,  and  source  for  convenient  reference  and  retrieval.  Parametric  fits 
were  performed  on  each  of  the  seven  terrain  types  to  reduce  the  data  to 
analytic  expressions.  These  results  are  discussed  in  Section  3.  In  general, 
the  fits  agree  qualitatively  with  those  performed  by  other  investigators. 

They  differ  from  others,  however,  in  that  they  cover  data  from  a variety 
of  sources  and  in  the  particular  formulation  of  the  models.  These  results 
are  summarized  in  Table  1. 

In  addition  to  models  for  average  a°,  the  variation  of  clutter  echos 
about  the  mean,  and  in  particular,  the  probability  of  them  exceeding  the  mean 
by  a great  deal,  is  of  great  concern  for  detection  of  targets  in  clutter. 

These  issues  are  discussed  in  Section  4.  There  the  different  mechanisms, 
variation  in  terrain  type,  local  undulations  in  slope,  man-made  objects, 
variation  in  local  conditions  (moisture  content,  orientation  of  crop  rows, 
etc.),  are  discussed  and  to  some  extent  quantified.  A theory  is  introduced 
for  predicting  the  pdf  of  clutter  based  on  the  underlying  pdf  of  the  indi- 
vidual scatterers,  resolution  cell  size,  and  spatial  autocorrelation  function. 
Based  on  the  theory,  examples  are  worked  out  for  five  terrain  types  and  a 
large  footprint  SBR  geometry  of  the  predicted  threshold  settings  for  two 
values  of  false  alarm  probability.  These  examples  show  that  the  thresholds 
differ  very  little  from  those  derived  from  a Rayleigh  model.  This  is  not 
to  say  that  clutter  statistics  are  always  Rayleigh  or  noise-like,  but  they 
should  be  much  more  so  for  large  resolution  cells  and/or  range  ambiguous 
geometries  than  for  otherwise.  This  is  the  first  time  scattering  theory 
(as  opposed  to  curve  fitting)  has  been  applied  to  statistics  of  SBR  clutter 
variation. 
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Clutter  spectra  and  their  impact  on  moving  target  indication  was 
discussed  in  Section  4 where  it  was  pointed  out  that  the  dominant  source 
of  clutter  spectral  broadening  was  due  to  platform  motion.  This  effect 
is  very  predictable  and  there  are  known  signal  processing  techniques  for 
illuminating  it  at  least  in  principle.  Internal  clutter  motion  is  usually 
smaller  but  harder  to  deal  with.  These  are  summarized  in  Table  7. 

Finally  in  Section  6 the  current  status  of  our  ability  to  model  clutter 
was  reviewed  and  shortcomings  in  the  present  data  base  were  pointed  out 
Given  these  areas  where  more  data  are  needed,  three  candidate  clutter 
measurement  programs  were,  analyzed  and  summarized  in  Table  8. 

In  final  summary,  this  report  contains  a model  for  clutter  as  it  is 
encountered  in  space-based  radar  systems.  It  is  based  on  the  best  experi- 
mental data  and  scattering  theory  currently  available.  Many  of  the 
approaches,  especially  in  Section  4 are  unique.  It  is  hoped  that  these 
results  will  be  useful  in  analyses  of  detection  of  targets  in  clutter  in 
general  as  a stepping-stone  for  other  modeling  efforts. 
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APPENDIX 

SIMULATION  OF  RADAR  SIGNALS 
A.l  CLUTTER  TAPE  GENERATION 


DSA  has  developed  a computer  model  which  simulates  the  time-dependent 
returns  from  ground  clutter  (in  digital  format)  as  seen  from  a space  radar 
system.  This  process  is  being  utilized  to  create  a data  tape  used  to 
evaluate  candidate  signal  processing  designs  for  RADC. 

The  user  must  specify  a number  of  parameters  and  option  choices.  The 
terrain  is  defined  by  a map  of  the  ground  clutter  as  observed  by  a particu- 
lar satellite  configuration,  also  specified.  The  radar  parameters  define 
antenna  pattern  and  waveform  including  beamwidth,  pulse  width,  sampling 
time,  and  number  of  pulses.  The  beam  is  divided  into  a specified  number  of 
ambiguous  regions.  Clearly,  the  numerous  variable  parameters  make  it  a 
versatile  program  with  modeling  capabilities  in  many  areas. 

The  block  diagram  in  Fig.  A-l  shows  basically  how  the  model  combines 
the  parameters  to  get  the  desired  signal  simulation.  The  backscatter  co- 
efficient map  will  include  the  different  terrain  types  in  the  scene  along 
with  the  model's  boundaries  and  inherent  spatial  correlations.  This  map 
along  with  the  radar's  antenna  pattern,  waveform,  and  target-sensor  geometry 
are  illustrated  in  Fig.  A-2.  These  parameters  in  combination  with  the 
orbital  kinematics  which  define  the  satellite's  altitude,  velocity,  and 
orientation,  are  input  to  clutter  integration  routines  to  obtain  a clutter 
range,  range-rate  map.  This  map  contains  signals  for  each  ambiguous  region 
folded  into  an  unambiguous  Doppler  region.  A model  for  the  random  fluctu- 
ation of  the  clutter  is  then  imposed  and  these  resultant  signals  are  used  to 
generate  a received  signal  for  each  range  strip.  The  signals  from  ambiguous 
range  strips  are  combined  and  then  convolved  with  a specified  single^ 
pulse  waveform.  Various  amounts  of  implementation  error  can  be  added  at 
this  point  for  parametric  analysis. 
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Figure  A-2.  Clutter  Map  Inputs 


A.  2 DOPPLER 


The  clutter  signal  as  a function  of  time  (t)  is  given  as 

C(t)  = cn  [ . A(R,R)e4lT,itR/X  s (t  - 2R/c)dRdR  (A-l) 

U-/R,R 

where  s(t)  is  the  waveform,  Cq  is  a complex  constant,  X is  the  wavelength, 

R is  the  range,  c is  the  speed  of  light,  R is  the  range  rate,  j = /^l,  and 
A{R,R)  is  the  complex  reflection  coefficient  of  clutter  whose  expected  value 

• p 

E|A(R,R)j  is  the  clutter  RCS  per  unit  range  and  range  rate*  The  latter, 
the  clutter  RCS  density,  or  "clutter  map"  is  a subject  in  itself  and  re- 
quires careful  modeling  to  faithfully  include  effects  of  viewing  geometry 
and  terrain  type.  It  requires  antenna  patterns,  platform  and  earth  motion, 
and  viewing  angles.  Equally  important  and  least  well-defined  is  the  sigma- 
zero  model  which  is  also  required.  Essentially  the  radar  cross-section  per 
unit  area  for  every  point  on  the  earth,  could  be  handled  by  maps  based  on 
actual  clutter  data,  a stochastic  model  based  on  probability  distributions 
and  spatial  correlations  of  actual  data,  or  a suitable  combination  of  the 
two.  Realistically,  the  last  solution  is  probably  the  most  suitable  in  most 
cases. 


For  notational  simplicity  we  define  the  following  quantities 


t*  = 2R/c 


v’  = 2R/X 


W 

W 


N 


e_2irj 

e-2*j/N 


Hereafter,  j can  be  used  as  an  index  as  well  as  since  its  meaning  should 
be  clear  from  context. 


For  the  moment  consider  the  waveform  to  be  a weighted  train  of  6 
function  pulses  - 
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N-l 


s(t)  - Sn  £ w 6(t  - nA) 
u n=0  n 


(A-2) 


where  Sg  is  a constant,  wn  is  a set  of  complex  weights  and  A is  the  inter- 
pulse period.  The  6-function  is  defined  such  that  for  an  arbitrary  func- 
tion g: 


J g(t)  6(t  - t* )dt  = g(t') 

— 00 

Substituting  (A-2)  and  (A-l)  and  integrating  yields 

oo 

c(t)  = CqSqS  / dv‘  A(t  - nA,  v')W"tv'wr 


(A-3) 


n -» 


Now  suppose  that  c(t)  is  sampled  at  discrete  times  tQ  + mA  for  N 
sample  (pulses). 


c(tg  + mA)  = cns 


030 


E/’d,- 


C0S0 


E/' 

n —co 


A(tg  + m - n)A,  v')W"mAv  wf 


dv'  A(tg  + nA,  v)W_rnAv  wm_n 


= 2 cn(t0  * m&)  - (N  - l)<n<N  - 1 


The  procedure  we  use  from  this  point  on  is: 


(A-4) 


1.  Generate  a(R,R),  based  on  a microscopic  description  of  a clutter 
scene,  i.e.,  you  specify  a(X,Y)  or  ^(latitude,  longitude)  and  use 
the  viewing  geometry  to  transform  it  into  o(R,R). 

2.  Apply  random  deviations  about  the  mean  to  obtain  c'(R,R)  where 

9 • 

the  expected  value  is  a(R,R).  Form  A(R,R)  is  such  that  the  phase 
is  random  and  o'  - A*A. 

3.  By  inverse  FFT  obtain  Cri(t)  for  a large  number  of  pulses.  [A 
pulse  is  a sample  of  the  discrete  series  Cn(ti).]  Truncate  to  an 
N pulse  burst. 
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4.  Use  Eq.  A-4  to  form  C(t0  + mA). 

We  repeat  for  a large  number  of  range  R values  and  convolve  with  the 
single-pulse  waveform  as  described  below. 

A.  3 PULSE  COMPRESSION  WAVEFORM 

This  procedure  so  far  is  valid  only  for  5-function  waveforms  such  as 
Eq.  A-2.  In  general  the  waveform  will  be  a train  of  identical  pulses  of 
finite  length.  This  is  given  by 

S(t)  = j~ dt"  s0(t")  T(t  - t") 

where  T(t)  is  a train  of  6-function  as  in  Eq.  A-3  and  sQ(t)  is  the  envelope 
of  the  individual  pulses. 

Equation  A-3  now  becomes 


c(t)  = c0s0 


jvf* 


A(t‘,  v')w 


Jdv'  £«(t  - r - t" 


- nA) 


X wns0^t' ' ^ 


C0S0  A(t  - t"  - nA*  v>rtv  wns0(f) 

- t‘ ' ).s0(t * ' )dt"  where  c^(t)  is  now  given 
by  Eq.  A-3. 


Since  the  process  of  convolving  a long  sequence 

cV  + m)  tm1n  < t°  < tmax,m  = 0,...,  N - 1 

with  a short  one  such  as  SqU)  (the  individual  pulse  length  is  certainly 
less  than  a and  probably  much  shorter)  is  straightfonvard,  it  seems  appro- 
priate to  first  compute  c^(t)  and  then  convolve  it  with  3Q(t)  by  means 
similar  to  the  FFT  convolution  algorithm. 

A.  4 RANGE-WALK 


The  effect  of  range  walk  is  a phenomenon  which  occurs  in  a high- 
resolution  radar  on  a fast  moving  platform.  We  restrict  this  discussion  to 
targets  whose  inherent  velocities  are  small  enough  so  as  not  to  allow 
additional  range-walk  in  relation  to  the  ground. 

The  general  expression  for  the  clutter  reflected  from  Earth  and  re- 
ceived by  the  radar  is  given  in  Eq.  A-l. 

* 

The  effect  of  range-walk  is  that  R itself  is  a function  of  R and  t. 
Im  fact  the  round  trip  delay  t equals 

t = 2R  c/c 
= 2RQt/c 

where  Rq  is  the  range  rate  of  targets  and  clutter  !n  the  center  of  the 
mainbeam.  The  approximation  obviously  does  not  hold  for  all  clutter,  how- 
ever we  are  primarily  interested  in  objects  in  the  mainbeam  only. 

Note  that  the  integral  over  R is  an  inverse  Fourier  transform. 
Therefore, define  for  a uniformly  spacial  (PRI  = A)  train  of  pulses 

r • , • 4TTjRt(,/A 

A'  (R,k)  = J dR'A(M)  e k (A-5) 
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where  0 < k <_  N - 1 , tk  3 kA  and  N is  the  number  of  pulses.  For  purposes 
of  this  discussion  we  ignore  such  effects  as  range  foldover  as  they 
are  not  relevant  to  this  discussion. 

Now  consider  the  return  from  the  Kth.  pulse.  We  make  the  additional 
approximation  that  the  range-walk  during  one  pulse  repetition  interval 
(PRI)  is  much  less  than  the  range  resolution. 


Rq  A « 6R  and  let  t = t + kA 


Ck(f)  = /dR  A,(R,k)S0(f  - R kA  - 2R/c) 


(A-6) 


From  here  on  we  define  A^( R)  s A'{R»k). 


Now,  the  fact  is  that  the  processor  is  digital  and  the  subpulse  wave- 
form is  sampled  at  discrete  times.  The  effect  is  that  corresponding  time 
samples  for  different  pulses  (different  k)  pertain  to  different  clutter. 


C.(t')  = J" dR  A'(R  ■*  R0kA,k)S0(f  - 2R/c) 


(A-7) 


In  the  simulation,  the  above  convolution  is  performed  discretely.  Define 
t"  = 2 R/C. 


Cfcttf,)  = 2 A-(t"  + t°)S0(t<  - t") 


(A-8) 


where  t?  = 2R  kA/c. 
k o 

Now  the  problem  is  that, in  generating. A(R,R);  it  is  hot  practi- 
cal to  divide  the  range  axis  into  a grid  much  finer  than  the  range 
resolution.  The  fact  that  the  change  in  tjj  from  one  pulse  to  the  next 
(2Rqa/c)  is  much  smaller  than  5t"  = 25R/c  where  6R  is  the  range  resolution. 
Therefore,  we  resort  to  an  interpolation. 


The  interpolation  can  be  difficult  because  A'  is  complex  and  can 
vary  greatly  from  one  sample  to  the  next.  After  considering  the 
alternatives  available  for  interpolation’,' it  was  decided  to  base  it  on 
the  Fourier  components  of  the  discrete  sample.  This  is  especially 
attractive  since  the  convolution  is  to  be  performed  by  the  convolution 
theorem  so  the  Fourier  components  are  required  anyway.  Let 


- > s1  t \ 1 1 


W'Z  K <V  e 


(A-9) 


Then  the  interpolation  gives 


1 

1 


? 
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As  we  have  said  this  memo  addresses  only  a procedure  for  range- 
walk  simulation,  not  compensation.  Obviously  the  perfect  way  to  com- 
pensate for  range  walk  in  this  simulated  signal  is  to  reverse  the 
process.  This  amounts  to  the  following 

1 VT  -2rrjf „ t?l  +2njf  t. 

Ak(ti)  = I ZLck(f£)  e £ kJS  (f^)  e * 1 (A-12) 

& 

-2ir jf  t . 

where  C^f*)  ■ 2.Ck(tf)  e 81 

This  again  amounts  to  interpolating  the  Ak(t)  so  that  the  same 
targets  and  clutter  appear  in  the  same  range  samples  for  every  pulse. 
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